Symposium on Theoretical Aspects of Computer Science 2009 (Freiburg), pp. 63-74 
www.stacs-conf.org 



o 
o 

(N 

D 



C/3 



(N 
> 

O 
00 

p 

O 
oo 

o 



X 



SHORTEST PATHS AVOIDING FORBIDDEN SUBPATHS 



MUSTAQ AHMED ^ AND ANNA LUBIW^ 

^ David R. Cheriton School of Computer Science, University of Waterloo, Canada 
E-mail address: m6ahmed@uwaterloo.ca 
E-mail address: alubiw@uwaterloo.ca 



Abstract. In this paper we study a variant of the shortest path problem in graphs: given 
a weighted graph G and vertices s and t, and given a set X of forbidden paths in G, find a 
shortest s-t path P such that no path in X is a subpath of P. Path P is allowed to repeat 
vertices and edges. We call each path in X an exception, and our desired path a shortest 
exception avoiding path. We formulate a new version of the problem where the algorithm 
has no a priori knowledge of X, and finds out about an exception x £ X only when a path 
containing x fails. This situation arises in computing shortest paths in optical networks. 
We give an algorithm that finds a shortest exception avoiding path in time polynomial in 
\G\ and |X|. The main idea is to run Dijkstra's algorithm incrementally after replicating 
vertices when an exception is discovered. 



1. Introduction 

One of the most fundamental combinatorial optimization problems is that of finding 
shortest paths in graphs. In this paper we study a variant of the shortest path problem: 
given a weighted graph G{V,E), and vertices s and t, and given a set X of forbidden 
paths in G, find a shortest s-t path P such that no path in X is a subpath of P. We 
call paths in X exceptions, and we call the desired path a shortest exception avoiding path. 
We allow an exception avoiding path to be non-simple, i.e., to repeat vertices and edges. 
In fact the problem becomes hard if the solution is restricted to simple paths [20]. This 
problem has been called the Shortest Path Problem with Forbidden Paths by Villeneuve and 
Desaulniers [22]. Unlike them, we assume no a priori knowledge of X. More precisely, we 
can identify a forbidden path only after failing in our attempt to follow that path. This 
variant of the problem has not been studied before. It models the computation of shortest 
paths in optical networks, described in more detail in the "Motivation" section below. Note 
that when we fail to follow a path because of a newly discovered exception, we are still 
interested in a shortest path from s to t as opposed to a detour from the failure point. This 
is what is required in optical networks, because intermediate nodes do not store packets, 
and hence s must resend any lost packet. 
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This paper presents two algorithms to compute shortest exception avoiding paths in 
the model where exceptions are not known a priori. The algorithms take respectively 
0{knlogn + km) and 0{{n + L)log(n + L) + m + dL) time to find shortest exception 
avoiding paths from s to all other vertices, where n = \V\, m = \E\, d is the largest degree 
of a vertex, k is the number of exceptions in X, and L is the total size of all exceptions. 

Our algorithm uses a vertex replication technique similar to the one used to handle non- 
simple paths in other shortest path problems [G, 22]. The idea is to handle a forbidden path 
by replicating its vertices and judiciously deleting edges so that one copy of the forbidden 
path is missing its last edge and the other copy is missing its first edge. The result is to 
exclude the forbidden path but allow all of its subpaths. The main challenge is that vertex 
replication can result in an exponential number of copies of any forbidden path that overlaps 
the current one. Villeneuve and Desaulniers [22] address this challenge by identifying and 
compressing the overlaps of forbidden paths, an approach that is impossible for us since we 
do not have access to X. Our new idea is to couple vertex replication with the "growth" 
of a shortest path tree. By preserving certain structure in the shortest path tree we prove 
that the extra copies of forbidden paths that are produced during vertex replication are 
immaterial. Our algorithm is easy to implement, yet the proof of correctness and the 
run-time analysis are non-trivial. 

1.1. Motivation 

Our research on shortest exception avoiding path was motivated by a problem in optical 
network routing from Nortel Networks. In an optical network when a ray of light of a 
particular wavelength tries to follow a path P consisting of a sequence of optical fibers, it 
may fail to reach the endpoint of P because of various transmission impairments such as 
attenuation, crosstalk, dispersion and non-linearities [12, 17]. This failure may happen even 
though the ray is able to follow any subpath P' of P. This non-transitive behavior occurs 
because those impairments depend on numerous physical parameters of the traversed path 
(e.g., length of the path, type of fiber, wavelength and type of laser used, location and 
gain of amplifiers, number of switching points, loss per switching point, etc.), and the effect 
of those parameters may be drastically different in P than in P' [2]. Forbidden subpaths 
provide a straight-forward model of this situation. 

We now turn to the issue of identifying forbidden paths. Because of the large number 
of physical parameters involved, and also because many of the parameters vary over the 
lifetime of the component [3] , it is not easy to model the feasibility of a path. Researchers 
at Nortel suggested a model whereby an algorithm identifies a potential path, and then 
this path is tried out on the actual network. In case of failure, further tests can be done 
to pinpoint a minimal forbidden subpath. Because such tests are expensive, a routing 
algorithm should try out as few paths as possible. In particular it is practically impossible to 
identify all forbidden paths ahead of time — we have an exponential number of possible paths 
to examine in the network. This justifies our assumption of having no a priori knowledge of 
the forbidden paths, and of identifying forbidden paths only by testing feasibility of a path. 

The shortest exception avoiding path problem may also have application in vehicle 
routing. Forbidden subpaths involving pairs of edges occur frequently ( "No left turn" ) and 
can occur dynamically due to rush hour constraints, lane closures, construction, etc. Longer 
forbidden subpaths are less common, but can arise, for example if heavy traffic makes it 
impossible to turn left soon after entering a multi-lane roadway from the right. If we are 
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routing a single vehicle it is more natural to find a detour from the point of failure when 
a forbidden path is discovered. This is different from our model of rerouting from s upon 
discovery of a forbidden path. However, in the situation when vehicles will be dispatched 
repeatedly, our model does apply. 

1.2. Preliminaries 

We are given an directed graph G(y, E) with n = \V\ vertices and m = \E\ edges where 
each edge e £ E has a positive weight denoting its length. We are also given a source vertex 
s € V, a destination vertex t € V, and a set X of paths in G. The graph G together with X 
models a communication network in which a packet cannot follow any path in X because 
of the physical constraints mentioned in Sec. 1.1. We assume that the algorithm can access 
the set X of forbidden paths only by performing queries to an oracle. Each query is a path 
P, and the oracle's response is either the confirmation that P is exception avoiding, or else 
an exception x G X that is a subpath of P and whose last vertex is earliest in P. Ties can 
be broken arbitrarily. In our discussion we say "we try a path" instead of saying "we query 
the oracle" because the former is more intuitive. In Sec. 4 we modify our algorithm for the 
case of an oracle that returns any exception on a path (not just the one that ends earliest). 
This requires more calls to the oracle but gives a faster run-time. 

We want to find a shortest path from s to i that does not contain any path in X as a 
subpath — we make the goal more precise as follows. A path is a sequence of vertices each 
joined by an edge to the next vertex in the sequence. Note that we allow a path to visit 
vertices and edges more than once. If a path does not visit any vertex more than once, we 
explicitly call it a simple path. A simple directed path from vertex v to vertex u; in G is called 
a forbidden path or an exception if a packet cannot follow the path from v to w because of 
the physical constraints. Given a set A of forbidden paths, a path ("^1,^2,^3, . . . ,vi) is said 
to avoid A if {vi,Vi+i, . . . Vj) A for all i,j such that I < i < j < I. A path P from s to i 
is called a shortest A-avoiding path if the length of P is the shortest among all A-avoiding 
paths from s to t. We will use the term "exception avoiding" instead of "X-avoiding" when 
A is equal to X, the set of all forbidden paths in G. 

1.3. Related work 

A shortest s-t path in a graph can be computed in 0{nlogn + m) time and linear space 
using Dijkstra's algorithm with Fibonacci heaps if all edge weights are non-negative, and 
in 0{mn) time and linear space using the Bellman- Ford algorithm otherwise [o]. When 
the edge weights are non- negative integers, the problem can be solved in deterministic 
0(771 log log n log log log n) time and linear space if the graph is directed [13] , and in optimal 
0{m) time if the graph is undirected [21]. In many of these cases, there are randomized 
algorithms with better expected times as well as approximation schemes. See Zwick [2.3] 
for a survey of shortest path algorithms, and Cabello [4], Goldberg and Harrelson [11] and 
Holzer et al. [15] for some of the more recent work. 

Two recent papers on shortest paths in graphs address the issue of avoiding a set 
of forbidden paths, assuming that all the forbidden paths are known a priori. The first 
paper gives a hardness result. Szeider [20] shows, using a reduction from 3-SAT, that the 
problem of finding a shortest simple exception avoiding path is NP-complete even when 
each forbidden path has two edges. If the forbidden paths are not known a priori, the 
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hardness result still applies to the case of simple paths because the lack of prior knowledge 
of the forbidden paths only makes the problem harder. 

The second paper, by Villeneuve and Desaulniers [22], gives an algorithm for a shortest 
(possibly non-simple) exception avoiding path for the case when all the forbidden paths 
are known a priori. They preprocess the graph in 0{{n + L) log(n + L) + m + dL) time 
and 0{n + m + dL) space so that a shortest path from s to a query vertex can be found 
in 0{n + L) time. They first build a deterministic finite automaton (DFA) from the set 
of forbidden paths using the idea of Aho and Corasick [1], which can detect in linear time 
whether a given path contains any of the forbidden paths. They then "insert" the DFA into 
G by replicating certain vertices of G in the manner introduced by Martins [(»], and then 
build a shortest path tree in this modified graph. Their algorithm cannot handle the case 
where the set of all forbidden paths is not explicitly given. Our algorithm is strictly more 
general, and we show in Sec. 4 that it solves their problem in roughly the same time but in 
less {0{n + m + L)) space. 

We now mention two problems that seem related to ours, but do not in fact provide 
solutions to ours. The first one is maintaining shortest paths in a dynamic graph, i.e., where 
nodes or edges may fail [7, 9, 14], or edge weights may change (e.g., [7, S]). Forbidden 
paths cannot be modeled by deleting edges or by modifying edge costs because all edges 
in a particular forbidden path may be essential — see Fig. 1 for an example. The second 
seemingly related problem is finding the k shortest paths in a graph. This was the subject 
of Martins [(i] who introduced the vertex replication technique that we use in our algorithm. 
There is considerable work on this problem, see Eppstein [10] for a brief survey. But the 
k shortest path problem is again different from our situation because a forbidden subpath 
may be a bottleneck that is present in all of the k shortest paths even for k G il(2"'^), see 
Villeneuve and Desaulniers [22]. 

In the context of optical networks researchers have studied many theoretical problems. 
See Ramaswami and Sivarajan [19] for details on optical networks, and Lee and Shay- 
man [17] and McGregor and Shepherd [18] for a brief survey of the theoretical problems 
that have been investigated. In the previous work, the effect of physical constraints on 
paths in optical networks is either not considered at all (e.g., Khuller et al. [16]), or simply 
modeled by a known constant upper bound on the length of such a path (e.g., Gouveia 
et al [12], Lee and Shayman [17] and McGregor and Shepherd [18]). To the best of our 
knowledge, none of the previous work on shortest paths in optical networks considers the 
fact that it is practically infeasible to know a priori all the forbidden paths in the network, 
i.e., all the constraints in X. Our paper handles the issue of physical constraints from a 
different and much more practical perspective. 

2. Algorithm for a shortest s-t path 

In our algorithm we begin with a shortest path tree rooted at s, ignoring the exceptions. 
We then "try out" the path from s to t in the tree. If the path is free of exceptions, we are 
done. Otherwise, to take the newly discovered exception into account, we modify the graph 
using path replication as described in the Introduction, and we modify the shortest path tree 
to match. In general, we maintain a modified graph and a shortest path tree in the graph 
that gives a shortest path in the original graph from s to every other vertex avoiding all the 
currently-known exceptions. We will first illustrate the idea with an example. Consider the 
graph G in Fig. 1(a), where the integers denote edge weights, and the dashed arrow marks 
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Figure 1: (a) Shortest paths and (b) shortest x-avoiding paths in a graph, where x 
{s,a,b,t). 



the forbidden path x = {s,a,b,t). Note that for simphcity we have used undirected edges 
in the figure to denote bidirectional edges. It is not hard to see that P = (s, c, a, b, t) is 
the shortest x-avoiding path from s to t. To find P, we first construct a shortest path tree 
rooted at s (marked using the heavy edges in Fig. 1(a)), and then try the path (s, a, b, t) in 
the tree. The path fails because it contains x, so we use a vertex replication technique similar 
to the one by Martins [G] to make duplicates of vertices a and b and delete edges (s, a') and 
(6, t), as shown in Fig. 1(b). We then construct a shortest path tree rooted at s (marked 
using the heavy edges in Fig. 1(b)) in the modified graph, and try the path {s,c,a',b',t) 
which "represents" the path P in G. We are done if x is the only forbidden path in G. 
Note that this approach can double the number of undiscovered forbidden paths. Suppose 
y = (c, a, b) is another forbidden path in G. We have two copies of y in the modified graph: 
{c,a,b) and {c,a',b'), and we have to avoid both of them. Our solution to this doubling 
problem is to "grow" the shortest path tree in such a way that at most one of these two 
copies is encountered in future. Our algorithm is as follows: 

1 construct the shortest path tree Tq rooted at s in Go = G; 

2 let i = 1; 

3 send a packet from s to t through the path in Tq; 

4 while the packet fails to reach t do 

5 let Xi be the exception that caused the failure; 

6 construct G, from Gj_i by replicating the intermediate vertices of Xi and then 
deleting selected edges; 

7 construct the shortest path tree Tj rooted at s in Gi using Tj_i; 

8 send a packet from s to t through the path in T^; 

9 let i = i + 1; 

In the above algorithm, the only lines that need further discussion are Lines 6 and 7; 
details are in Sections 2.1 and 2.2 respectively. In the rest of the paper, whenever we focus 
on a particular iteration i > 0, we use the following notation: (i) the path from s to t in 
Tj_i, i.e., the path along which we try to send the packet to t in Line 4 in the iteration, is 
(s, vi,V2, ■ ■ ■ , Vp, t), and (ii) the exception that prevented the packet from reaching t in the 
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(b) G;_i 



(c) Gi 



Figure 2: Modifying Gj_i to Gj: (a) The part of Gj_i at an exception 
(vr-3,Vr-2,Vr-i,Vr,Vr+i), with / = 3. (b) Rephcating vertices to create 
G^_^. The dashed paths show two of the 8 copies of the exception, (c) Deleting 
edges to create Gj. The dotted hnes denote deleted edges. 

2.1. Modifying the graph 

The modification of Gj_i into Gi (Line 6) in the ith iteration eliminates exception Xi 
while preserving all the Xj-avoiding paths in Gi-i. We do the modification in two steps. 

In the first step, we create a graph G^_i by replicating the intermediate vertices of Xi 
(i.e., the vertices Vr-i+i,Vr-i+2, ■ ■ ■ ,'Ur)- We also add appropriate edges to the replica v' 
of a vertex v. Specifically, when we add v' to Gj_i, we also add the edges of appropriate 
weights between v' and the neighbors of v. It is easy to see that if a path in Gj_i uses /' < I 
intermediate vertices of Xi, then there are exactly 2 copies of the path in G[_i. We say 
that a path in G^_^ is Xj-avoiding if it contains none of the 2 copies of Xj. 

In the second step, we build a spanning subgraph Gi of G'^_i by deleting a few edges 
from G^_^ in such a way that all copies of Xi in G[_i are eliminated, but all Xj-avoiding 
paths in G[_i remain unchanged. To build Gi from G[_i, we delete the edges (fj_i,u') 
and {v'-,Vj^i) for all j € [r — I + l,r]. We also delete the edge {vr,Vr+i), all the outgoing 
edges from -y^ except (v'^^Vr+i), and all the outgoing edges from f' except (f',f'_^;^) for all 
j S [r — / + 1, r — 1]. Figure 2 shows how the "neighborhood" of an exception changes from 
Gj_i to Gi. As before, the undirected edges in the figure are bidirectional. 

Observation 2.1. Graph Gi has no copy of Xj. 

In Sec. 3.1 we will prove that Gi still contains all the Xj-avoiding paths of Gj_i. 

The vertices in Gi [G[_-^] that exist also in Gj_i (i.e., the ones that are not replica 
vertices) are called the old vertices of Gi [respectively G[_i]. Note that the vertices of Go 
exist in Gj for all i > 0. These vertices are called the original vertices of Gi. 



2.2. Constructing the tree 

In Line 7 of our algorithm we construct a tree Tj that contains a shortest Xj-avoiding 
path from s to every other vertex in Gj_i. Tree Tj is rooted at s, and its edges are directed 
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away from s. Not every shortest path tree rooted at s in Gi will work. In order to guarantee 
termination of the algorithm, Tj must be similar to Tj_i, specifically, every Xj-avoiding path 
from s in Tj_i must be present in Tj. The necessity of this restriction is explained in Sec. 3.3. 
We construct the required Tj by preserving as much of Tj_i as possible. We apply 
Dijkstra's algorithm starting from the part of Ti_i that can be preserved. Let V be the 
set of vertices that are either replica vertices in Gi, or descendants of Vr+i in Tj_i. We first 
set the weight of each f S y to infinity, and temporarily set Tj = Tj-i — V'. Then, for 
each V G V, we set the weight of v to the minimum, over all edges (u, v), of the sum of the 
weight of u and the length of (u, v). Finally, we initialize the queue used in Dijkstra with 
all the vertices in V' and run the main loop of Dijkstra's algorithm. Each iteration of the 
loop adds one vertex in V' to the temporary Tj. When the queue becomes empty, we get 
the final tree Tj. 

3. Correctness and analysis 

3.1. Justifying the graph modification 

In this section we prove the following lemma, which uses the notion of a corresponding 
path. Consider any path Pi in Gi. By substituting every vertex in Pi that is not present 
in Gi_i with the corresponding old vertex in Gj_i, we get the corresponding path Pi^i in 
Gj-i. This is possible because any "new" edge in Gi is a replica of an edge in Gj^i. We 
define the corresponding path Pj in Gj for all j < i hy repeating this argument. 

Lemma. If Pi is a shortest path from s to an original vertex v in Gi, Pq is a shortest 
{xi, X2, . . . , Xj }-avoiding path from s to v in Gq. 

To prove the above lemma (repeated as Lemma 3.3 below), we will first prove that Xj- 
avoiding paths in Gj_i are preserved in Gi (Lemma 3.2), using the following characteristic 
of an Xj-avoiding path in the intermediate graph G^^^: 

Lemma 3.1. For any Xi-avoiding path P from s to v that uses only the old vertices in G^_i, 
there exists a copy of P in Gi that starts and ends at the old vertices s and v respectively, 
and possibly passes through the corresponding replicas of its intermediate vertices. 

Proof. Graph Gi contains all the edges between pairs of old vertices in G^_^ except for the 
directed edge {vr,Vr+i). Thus P can remain unchanged if it does not use this directed edge. 
Otherwise we will re-route any portion of P that uses the directed edge {vr,Vr+i) to use 
the replica edge {v'j.,Vr+i) instead. Let P = {s = wi,W2, ■ ■ ■ ,Wq = v), and {wj,Wj+i) be 
an occurrence of {vr,Vr+i) in P. Tracing P backwards from Wj, let h < j he the minimum 
index such that {wh, Wh+i, . . . , Wj+i) is a subpath of Xj. Because P is Xj-avoiding, Wh must 
be an intermediate vertex of Xj . This implies that h > 1, since s = wi is not an intermediate 
vertex of Xj because of the following reasons: (i) Xj is a path in the shortest path tree rooted 
at s in Gj, and (ii) there is no replica of s in Gj. Therefore Wh-i exists. We will reroute the 
portion of P between w^^i and Wj+i by using the corresponding replica vertices in place of 
the subpath {wh , . . . ,Wj) of Xj . Note that the required edges exist in Gj (since P does not 
contain the whole exception Xj), and that the portions of P that we re-route are disjoint 
along P. Moreover, P starts and ends at the old vertices s and v respectively. ■ 
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Lemma 3.2. Any Xi-avoiding path from s to v in Gj-i has a copy in Gi that starts and 
ends at the old vertices s and v respectively, and possibly goes through the corresponding 
replicas of its intermediate vertices. 

Proof. Let P be the Xj-avoiding path in Gj-i. As we do not delete any edge to construct 
G^_i from Gi-i, P remains unchanged in G[_^. Moreover, P uses no rephca vertex in G[_]^. 
So, Lemma 3.1 imphes that P exists in Gi with the same old vertices at the endpoints, 
possibly going through the corresponding replicas of the intermediate vertices. ■ 

Lemma 3.3. If Pi is a shortest path from s to an original vertex v in Gi, Pq is a shortest 
{xi, X2, • • • , Xi} -avoiding path from s to v in Gq. 

Proof. For any j G [0,i], let Xj = {xj+i,Xj+2; • • • ,Xi}. We show that for any j, if Pj is a 
shortest Xj-avoiding path in Gj, then Pj-i is a shortest Xj_i-avoiding path in Gj-i. The 
lemma then follows by induction on j, with basis j = i, because Xj = and thus Pi is a 
shortest Xj-avoiding path in Gi. 

If Pj is a shortest Xj-avoiding path in Gj, Pj is Xj_i-avoiding because Pj is Xj-avoiding 
by Observation 2.1, and Xj U {xj} = Xj-i. So, the corresponding path Pj-i is also Xj^i- 
avoiding. If we assume by contradiction that Pj~i is not a shortest Xj_i-avoiding path in 
Gj^i, then there exists another path P'^_i from s to f in Gj-i which is Xj_i-avoiding and 
is shorter than Pj-i- Since Xj G Xj-i, P'j_i is Xj-avoiding, and hence by Lemma 3.2, there 
is a copy P' of path Pj_i in Gj which has the same original vertices at the endpoints. As 
Pj^i is Xj-avoiding, Pj is also Xj-avoiding. This is impossible because P' is shorter than 
Pj. Therefore, Pj-i is a shortest Xj_i-avoiding path in Gj^i. m 

3.2. Justifying the tree construction 

To show that the "incremental" approach used in Sec. 2.2 to construct Tj is correct, we 
first show that the part of Tj„i that we keep unchanged in Tj is composed of shortest paths 
in Gi-. 

Lemma 3.4. For every vertex v that is not a descendant of Vr+i in Tj_i, the path P from 
s to V in Tj-i is a shortest path in Gi. 

Proof. First we show that P exists in Gi. Every vertex in Tj-i exists in Gi as an old vertex. 
So, P exists in Gj through the old vertices if no edge of P gets deleted in Gj. The only 
edge between a pair of old vertices in Gj_i that gets deleted in Gj is {vr,Vr+i). Since v is 
not a descendant of Vr+i in Tj_i, P does not use the edge {vr,Vr+i). Therefore, no edge of 
P gets deleted in Gj. So, P exists in Gj through the old vertices. 

Neither the modification from Gi_i to G^_]^ nor the one from G[_-^ to Gj creates any 
"shortcut" between any pair of vertices. So, there is no way that the distance between a 
pair of old vertices decreases after these modifications. Since these modifications do not 
change P, which is a shortest path in Gj_i, P is a shortest path in Gj. ■ 

Lemma 3.5. The tree Ti is a shortest path tree in Gi. 

Proof. For every vertex v that is not a descendant of Vr+i in Tj-i, the path P from s to v 
in Ti is the same as the one in Tj-i and hence, a shortest path in Gj (Lemma 3.4). For all 
other vertices u in Gj, it follows from Dijkstra's algorithm that the path from s to « in Tj 
is a shortest path. ■ 
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Lemmas 3.3 and 3.5 together prove that our algorithm is correct provided it terminates, 
which we estabhsh in the next section. 

3.3. Analyzing time and space requirement 

Although in every iteration we eliminate one exception by modifying the graph, we 
introduce copies of certain other exceptions through vertex replication. Still our algorithm 
does not iterate indefinitely because, as we will show in this section, the incremental con- 
struction of the shortest path tree (Sec. 2.2) guarantees that we do not discover more than 
one copy of any exception. We first show that any exception in Gi_i has at most two copies 
in Gi (Lemma 3.6), and then prove that one of these two copies is never discovered in the 
future (Lemma 3.7): 

Lemma 3.6. Let y ^ x, he any exception in Gi^i. If the last vertex of y is not an 
intermediate vertex of Xi, then Gi contains exactly one copy of y. Otherwise, Gi contains 
exactly two copies of y. In the latter case, one copy of y in Gi ends at the old vertex v and 
the other copy ends at the corresponding replica v' . 

Proof. Let vr = {wi,W2, . . . ,Wj) he a maximal sequence of vertices in y that is a subsequence 
of {vr~i-\^i,Vr-i+2, ■ ■ ■ , ^r)- Let w'j be the replica of Wj in Gi. We will first show that if there 
is a vertex f in y right after vr, then exactly one of the edges (wj, v) and {w'j,v) exists in Gi. 
Consider the subgraph of Gi induced on the set of replica vertices {f^_;_,_i, v^_;^2' ■ ■ ■ ,K}'- 
this subgraph is a directed path from w^_;_|_i to v',., and the only edge that goes out of 
this subgraph is {v'j,,Vr+i). Therefore, (i) when {wj,v) = {vr,Vr+i), {w'j,v) G Gi and 
{wj,v) ^ Gi, and (ii) otherwise, {wj,v) € Gi and {w'-,v) Gi. 

Now Gi has exactly two copies of vr: one through the old vertices, and another through 
the replicas. The above claim implies that when there is a vertex v in y right after tt, Gi 
has at most one copy of the part of y from wi to v. However, when vr is a suffix of y, Gi 
has both the copies of the part of y from wi to Wj. The lemma then follows because any 
part of y that contains no intermediate vertex of Xi has exactly one copy in Gi. m 

Lemma 3.7. Let y ^ Xi he any exception in Gi^i such that the last vertex of y is an 
intermediate vertex v of Xi. The copy of y that ends at the old vertex v in Gi is not 
discovered hy the algorithm in any future iteration. 

Proof. The copy of the path (s, vi,V2, . . . , Vj-) through the old vertices in Gi contains v. Let 
P be the part of this path from s to v. Clearly, P € Tj_i, and P does not contain any 
exception because the oracle returns the exception with the earlier last vertex. So, the way 
we construct Tj from Tj_i for any iteration j >i ensures that P GTj. 

Let yi be the copy of y that ends at v. Now yi is not a subpath of P because P does 
not contain any exception. For any j > i, P € Tj, and both P and yi end at the same 
vertex, therefore yi ^ Tj. So, a packet in iteration j will not follow yi, and yi will not be 
discovered in that iteration. ■ 

Lemma 3.8. The while loop iterates at most k = \X\ times. 

Proof. For any iteration i, Gj_i contains Xj, and Gi does not contain Xj. Every exception 
other than Xj in Gj_i has either one or two copies in Gi (Lemma 3.6). By Lemma 3.7, if an 
exception has two copies in Gi, only one of them is relevant in the future. Thus the number 
of exceptions effectively decreases by one in each iteration. The lemma then follows. ■ 
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To determine the running time, observe that the number of vertices increases in each 
iteration. However, we run Dijkstra's algorithm on at most n vertices in any iteration, be- 
cause the number of rephca vertices added in each iteration is always less than the number 
of vertices in the part of the shortest path tree that is carried over from the previous tree 
in our incremental use of Dijkstra. Moreover, we can make sure that Dijkstra's algorithm 
examines at most m edges in iteration i, by deleting a few more edges from Gi after per- 
forming the graph modification described in Sec. 2.1. More precisely, for each old vertex 
V G {vr-i+i, Vr-i+2, • • • > ^r}, since the label (i.e., the "distance" from s) put on v by Dijk- 
stra's algorithm in the previous iterations remains unchanged later on, we can safely delete 
from Gi all the incoming edges of v without affecting future modifications. (Note that for 
all j € [r — l + l,r], old vertices Vj and Vj+i are no longer adjacent in Gi, although the edge 
{vj,Vj+i) still exists in Tj.) It is not hard to see that the number of new edges in Gi is now 
equal to the number of edges deleted from Gj_i. 

Theorem 3.9. The algorithm computes a shortest X -avoiding path in 0{knlogn + km) 
time and 0{n + m-\- L) space. 

Proof. The correctness of the algorithm follows from Lemmas 3.3 and 3.5. 

Let li be the number of intermediate vertices of the exception discovered at the ith. 
iteration (thus the size of the exception is li + 2). The ith iteration adds li vertices. Since 
the algorithm iterates k times (Lemma 3.8), there are n + ^j=^ li < n + L vertices in the 
graph at termination. Because in each iteration the number of added edges is equal to the 
number of deleted edges, the space requirement is 0{n + m + L). 

Each iteration of our algorithm takes 0(|^| log \V\ + \E\) = 0(n logn -|- m) time, and 
the total time requirement follows. ■ 

We note that in practice, the algorithm will not discover all k of the forbidden paths. 
It will discover only the ones that "interfere" in getting from s to t. 

4. Extensions 

This section contains: (1) an algorithm to compute shortest paths from s to every other 
vertex in G; (2) an analysis in the case when X is given explicitly; and (3) a version of 
the algorithm where the oracle returns any exception on a query path, rather than the 
exception that ends earliest. 

The algorithm in Sec. 2 can be extended easily to compute a shortest path from s to 
every other vertex in G. We simply repeat the previous algorithm for every vertex in G, but 
with a small change: in every iteration (except of course the first one) we use the graph and 
the shortest path tree constructed at the end of previous iteration. Since every exception 
in X is handled at most once, the Avhile loop still iterates at most k times, and therefore, 
the time and space requirements remain the same. 

Theorem 4.1. The algorithm computes shortest X -avoiding paths from s to all other ver- 
tices in 0{knlogn -\- km) time and 0{n + m + L) space. 

Our algorithm applies when X is known explicitly; taking into account the cost of 
sorting X so that we can efficiently query whether a path contains an exception we obtain: 

Theorem 4.2. When X is known a priori, we can preprocess the graph in 0{knlog{kn) -\- 
km) time and 0{n + m + L) space so that we can find a shortest X -avoiding path from s 
to any vertex in 0{n + L) time. 
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Recall that Villeneuve and Desaulniers [22] solved this problem in 0((n + L)log(n + 
L) + m + dL) preprocessing time, 0{n + m + dL) space and 0{n + L) query time. Our 
algorithm is more space efficient than theirs. Our preprocessing is slightly slower in general, 
although it is slightly faster in the special case L = @(kn) and m = o{dn) (intuitively, when 
the exceptions are long, and the average degree of a vertex is much smaller than the largest 
degree) . 

Finally, returning to the case where X is not known a priori, we consider a weaker 
oracle that returns any exception on the query path, rather than the exception that ends 
earliest. At the cost of querying the oracle more often, we obtain a better run-time. The 
idea is to query the oracle during the construction of a shortest path tree. The algorithm is 
very similar to Dijkstra's, the only difference is that it handles exceptions inside Dijkstra's 
loop. More precisely, right after a vertex v is dequeued and added to the current tree, we 
try the s-v path in the tree. If the path is exception avoiding, we update the distances of 
the neighbors of v and go to the next iteration, as in "traditional" Dijkstra's algorithm. 
Otherwise, we remove v from the current tree, perform vertex replication and edge deletion 
as described in Sec. 2.1, and then go to the next iteration. 

Theorem 4.3. The algorithm described above computes shortest X-avoiding paths from, s 
to all other vertices in 0{{n + L) log(n + L) + m + dL) time and 0{n + m + L) space. 

Proof. There are at most n + L vertices in the modified graph in any iteration. So, the loop 
in the modified Dijkstra's algorithm executes at most n + L times, and the priority queue 
holds at most n + L entries. Moreover, within Dijkstra's loop vertex replication and edge 
deletion take 0{dL) time in total. The running time then follows. The proof of correctness 
is similar to that of Theorem 4.1 except that the "current" shortest path tree is no longer 
a spanning tree in the current graph. ■ 

This new algorithm is faster than the old algorithm of Theorem 4.1 in general but makes 
as many as n + L queries to the oracle versus at most k oracle queries for the old algorithm. 
The old algorithm is slightly faster in the special case L = Q{kn) and m = o{dn). 

5. Conclusion 

Motivated by the practical problem of finding shortest paths in optical networks, we in- 
troduced a novel version of the shortest path problem where we must avoid forbidden paths, 
but we only discover the forbidden paths by trying them. We gave an easily implement able, 
polynomial time algorithm that uses vertex replication and incremental Dijkstra. 

As we have mentioned before, in practice our algorithms will not discover all the forbid- 
den paths in X. In fact, the running time of each of our algorithms is determined by only 
the forbidden paths that "interfere" in getting from s to t. An interesting open problem 
is to bound the number of such paths. We conjecture that in a real optical network, the 
number of such paths is o{k), and therefore, our algorithms run much faster in practice. 
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